YerPhI - 1554(3) - 2000 
NRCPS-HE-2000-21 



Direct CP-asymmetry in Inclusive Rare B-decays in 
o: 2HDM. 

O 

Q 

H. H. Asatryan, H. M. Asatriaif\, G. K. Yeghiyan 

> ; 

\jT) , Yerevan Physics Institute, 



2 Alichanian Br. St. Yerevan, Armenia 
G.K.Savvidy 



00 
O 
(N 

o 
o 

National Research Center Demokritos, 
Qh| Agia Paraskevi, GR-15310 Athens, Greece 

43 

> 

X 
S3 

Abstract 

The direct CP-asymmetry in the inclusive B — > AVy and B — > Xrfe + e~ decays is investi- 
gated in the two-Higgs doublet extension of the Standard Model (2HDM). The investigation 
is performed in the lowest non-vanishing order of the perturbation theory using the existing 
restrictions on the 2HDM parameters space. It is shown that the direct CP-asymmetry in the 
B — > Xffi decay can deviate significantly from the Standard Model predictions. In the presence 
of only one source of CP- violation (the CKM matrix weak phase) acp{B — > Xffi) can have the 
sign opposite to that in the SM. The new source of CP-violation can make \acp{B — > X^7)| 
arbitrary small (unlike the SM case) and hence unmeasurable. Quantitatively, the obtained re- 
sults suffer from the uncertainty of the choice of renormalization scale. As for the B — > Xae + e~ 
rate asymmetry, its renormalization scale dependence in the lowest non-vanishing order does 
not allow to conclude if this quantity is efficient for testing New Physics beyond the Standard 
Model. 
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1 Introduction 



The experimental evidence for the Standard Model Higgs boson |IJ increases the interest to the 
Higgs sector of the electroweak theory and in particular to the Two-Higgs doublet extension of 
the Standard Model. The study of the rare B-decays is known to be an important source of 
information for such an extension of the Standard Model and on the physics at the scales up to 
several TeV. In particular these decays combined with the studies of other CP-violating effects 
in K- and B-physics can be useful for understanding the CP-breaking phenomenon. 

The first experimental evidence for the exclusive rare decay channel B — > K* / y has been 
reported by the CLEO collaboration 0. It was followed by the observation of the inclusive 
decay B — > X s 7 B. At the present the experimental prediction for the B X s 7 branching 
ratio from CLEO is 

< BR{B -> X sl ) >= - (BR{B -> X sl ) + BR(B -> X sl )) = (3.15 ± 0.54) x 10~ 4 (1.1) 

For the other rare B-decays, such as B — > X^, B — > X s ^£ + £~ etc., only upper experimental 
bounds exist. Hopefully these decays will be observed at the B-factories which have already 
started operating. 

Huge theoretical efforts have been devoted to the investigation of the rare inclusive decays 
in the Standard Model (SM) and its extensions. The most thorough study has been carried 
out for the B — > X s 7 decay []5|-[28|. In particular, in [I5||-p2|| the NLO corrections have been 



calculated in the SM reducing the uncertainty in theoretical prediction for BR(B — > X s "f) to the 
level of a few percent. In the extensions of the SM the NLO corrections have been considered in 



Z3"|-||27||. In these models the theoretical uncertainty due to renormalization scale dependence 
is still significant. The theoretical prediction for BR(B — > X s ^) in the SM is in reasonable 
agreement with ( |1 . 1| ) . Consequently strong constraints on the parameters of the SM extensions 
can be derived. In particular such constraints have been obtained for the Two-Higgs Doublet 



Model (2HDM) ||, gj and the Minimal Supersymmetric Standard Model (MSSM) |5|-|7|. 

Another popular rare B-decay mode, B — > X s £ + £~, has been considered in the context of 
the SM and its extensions in || ||, [29|-[32|. It is expected that this decay mode will be also 



useful for testing the New Physics. 

The direct CP-asymmetry in rare B-decays is another important source of information on 
the physics beyond the SM. In the SM the only source of CP-violation is the single phase 
in the CKM-matrix. Due to known hierarchy of the CKM-elements, b — > s transitions have 
larger branching fractions while in b — > d transitions large CP-violating effects are expected. 
Numerically (in the SM) the branching ratio for the B — ► Xj'j is 10 -v- 20 times smaller than 
for B — > Xgj but the CP-asymmetry of the latter can reach 27% to be compared with the 
prediction —(0.4 4- 1)% |33j for the former. As it was pointed out in P4]| , the B — ► X^ rate 
asymmetry is statistically more accessible than that for B A s 7. A similar hierarchy between 
the branching ratios and asymmetries holds (in the SM) also for the decays B — > X s e + e~ and 



B^X d e + e- 32 



The aim of this paper is to perform a detailed study of the direct CP-asymmetry in the 
inclusive B — > X^7 and B — > X d e + e~ decays in the lowest non- vanishing order of the per- 
turbation theory within the general 2HDM - the simplest extension of the SM where another 
source of CP- violation is present. The b —>■ d'j transition asymmetry in the 2HDM has been 
examined in fl35| . In the present paper we improve the results of ref. fl35| , taking into account 



new theoretical results and experimental data. In particular the scale parameter dependence 
is investigated in details. The included effects lead to a significant changes in the predictions 
for the B — > X^ decay rate asymmetry: it reaches -25% when being negative and 37% when 
being positive. The CP-asymmetry in B — ► Xde + e~ decay in the SM and 2HDM has been 
studied in []3I| [J^] and |R7| respectively. In ref. |37| it was obtained that for a specific choice of 



the parameters of the theory, B — ► X d e + e~ decay asymmetry is more efficient for testing the 
physics beyond the SM than the branching ratio. In present paper we show that when taking 
into account the renormalization scale dependence in the lowest non-vanishing order the allowed 
interval of the CP-asymmetry for B — > X d e + e~ decay in the SM and 2HDM does not differ 
significantly. Therefore to clarify whether the CP-asymmetry can be useful for testing the va- 
lidity of the Standard Model one should include higher-order corrections to a C p(B — > X d e + e~) 
which have not been calculated yet. 

The paper is organized as follows. In Section 2 we briefly describe the 2HDM, give the 
formulae which are needed for the calculation of B — > Xjpf and B — > Xde + e~ rate asymmetries 
and discuss the existing restrictions on the parameters of the model. In Sections 3 and 4 we 
present the numerical results for acp(B — > X^) and a C p(B — > X d e + e~) respectively. Finally 
Section 5 is devoted to the summary of the obtained results and conclusions. 



2 Formulae for B —> X^7, B — > X^e + e Rate Asymmetries 
and Constraints on Parameters of the Model. 

2.1 The Two-Higgs Doublet Model. 

It is known that in the SM the rare decays B — > X d ( s )l and B — > X d ( s )e + e~a,re loop-induced: in 
the lowest order they proceed via exchange of up-type quarks and W ± bosons in the loops. In 
the 2HDM there are additional diagrams with W-bosons replaced by the charged Higgs bosons 
{H^). The Yukawa interaction of the quarks with the Higgs doublets $i and $2 m general is 

^ £ = g L (r?$i + r d 2 $ 2 )d R + q l {t^ x + r%$ 2 )u R , (2.1.1) 

where Ql is the left-handed quark doublet, ur and dp are the right handed quark singlets, 
$1,2 = 202^1,2 (02 is the Pauli matrix) and Tf, rf, r", T% are matrices in the flavor space. One 
usually considers two versions of the 2HDM: 

= Tf = (Model I); 

r^ = r^ = o (Model ii). 

In this paper we also consider (neglecting FCNC's of the Lagrangian) the general case of the 
2HDM (Model III), where all the quantities Tf, Tf, T", are nonzero, thus investigating 
possible impact of the Higgs doublets vacuum phase on the B — ► X^ and B — > X d e + e~ rate 
asymmetries. Note also that the results obtained in the Model III in a large extent are valid 
for the multi-Higgs doublet models with only one light charged Higgs particle |24] . 



The physical charged Higgs state is related to the charged components of $1 and $2 as |E 

$± = e ±i5 H ± sin /3, ®f = H ± cos (3 



■1 



where tan/? = v 2 /vi is the ratio of vev's of the Higgs doublets and 5 is the Higgs doublets 
vacuum phase which yields spontaneous CP-violation. The couplings of the charged Higgs 
boson to t R and b R have the following form: 

(v* d \ ( V, 

£ t m t H-{d, s, b) L V* | t R = ve i5 *{h tl sin/ftr 1 * - h t2 cos P)H~(d, s, b) L 

V v t t 
( v td 



£ b m b H + (u, c,t) L 



V ts | b R = ve~ i5b {h bl sin[3e iS - h b2 cos (3)H+(u, c,t) L 

V v tb 




V ts I b R (2.1.2) 

V v tb 



where ht 



(^2)33, h bl2 = (rf )2 ) 3 3, v = + v\ = 174GeV, V is the CKM matrix, m t and 
m b are the running top and bottom masses (for the relations between the quarks running and 

and 5 t , 5 b are related to the parameters of the theory by J5S 



pole masses see e.g. 



m t e 



m b e 



(h tl cos (3e %s + h t2 sin f3)v 
= (h bl cos f3e lS + hb 2 sin j3)v 



We see that in the most general version of the 2HDM the couplings £t and are complex so that 
there is a new source of CP-violation connected with spontaneous CP-breakdown. However in 
the Models I and II these parameters are real: £ t = = — cot (3 and £ t = — cot /3, = tan /3 
respectively. 

Some restrictions on the couplings h tl , h t2 , h bl , h b2 can be derived [O] assuming the absence 
of new physics between the electroweak breaking scale fj,w ~ lOOGeV and unification scales 
fiG ~ (10 15 — 10 19 )GeV and requiring the validity of the perturbation theory at the whole energy 
range between fi w and [Iq. As a result one obtains that at the electroweak breaking scale h tl , 
ht 2 , h bl , h b2 are of the order of unity or smaller. For the present experimental values of the 
top mass, mF t ole = (174 ± 5)GeV it was found that tan/? > 1 for the Models I and II. It 
is easy to see from (2.1.2) that |£t| ~ 1 or smaller. On the contrary one finds that |£&| can be 
much larger than unity when tan/3 ^> 1 or (in the Model III) tan/3 1. Strong restrictions on 
£t and can also be derived from the experimental constraints on BK(B — > X s 7) and B — B 
mixing. We discuss the experimental bounds on the parameters of the theory in subsection 2.3. 



2.2 Formulae for B — > Xd'J, B — > X^e + e decay asymmetries. 

We study the decays B — > Xm s )J, B — > X d ^e + e~ using the effective theory with five quarks 
obtained by integrating out the heavy degrees of freedom which are W and Z bosons, t-quark 
and the charged Higgs boson. The effective Hamiltonian for the decays B — > X^s)! can be 
written as 



H eff (b^d(sh(+g) 



4G F 
1/2 



ECi(M)0 4 (/O - A?'> £Ci(/i)(G?(/*) - 0,(a*) 

i=l »=1 



(2.2.1) 



2 If the 2HDM originates as an "effective" low-energy theory, then similar restrictions on the coupling con- 
stants of the theory can be obtained, assuming the absence of new physics between the "full" theory scale and 
unification scales. 



1 



where X^ = V t * d ^V t b, X^ = V^r^Vyj,. Here we present only the most important operators 

Ox = d(s) L YT a c L c Ll ,T a b L , lu = d(s) L YT a u L u Ll ,T a b L , 

2 = d(s) L y»c L c L ^b L , 2u = d{s) L l^u L u L y^b L) (2.2.2) 



7 = j^m b (f,) d(s) L a^b R F, v , 8 = j^Mti d(s) L a^T a b R G%. 



The remaining operators contribute mainly through the mixing and otherwise can be safely 
neglected. The full operator basis together with the expressions for the corresponding Wilson 



coefficients Ci can be found elsewhere [22, 24, 28 



The processes B — > Xd( s )l can be well described by the partonic level transitions b — > d(s)y 
and b — > d{s) r yg. The calculation of these partonic transitions includes the following three steps 
TBI: 



1. The Wilson coefficients Ci at the heavy mass scale fiw ~ Mw must be calculated, match- 
ing the effective and full theories. In the next-to-leading order the matching has to be 
done at the 0(a s ) level, i.e. C^w) = Cf\nw) + a sl (^)C\ 1 \ijlw)- 

2. The renormalization group equations (RGE) must be used to obtain the Wilson coeffi- 
cients at the low-energy scale /i ~ m^. In the next-to-leading order this step requires the 
knowledge of the anomalous dimension matrix up to the order a 2 s . 

3. The matrix elements of the operators Oi for the processes b — * d(s) r y, b — > d(s)yg have to 
be calculated. 

Since the effective Hamiltonian (2.2.1) is the same for the 2HDM and SM, the 2HDM effects 
enter only through the first step. At the 0(a s ) level the matrix elements of Oi for the processes 



b — > 57 and b — > syg have been calculated in ref's [|T^, 16, 17, 21]. The obtained results are 
relevant also for b — > dy and b — > dyg decays when making an obvious replacement A| — > Xf , 
i=t,u. At the order a 2 the anomalous dimension matrix can be found in m| . The explicit 
relations between the Wilson coefficients at low and high energy scales are given in |Z2], |24|| . 
Here are the leading order expressions for the Wilson coefficients, corresponding to the subset 
of the relevant operators (2.2.2) 

C?V) = ^/23_ r 12/23 ) 

cfV) = *f 6/23 ^% w ) + ^(^ (2.2.3) 



where fj = a s (//w)/a s (/x), and the numerical values of hj, dj, hp a - are available in p2| ^4 . 

It is known that in the leading logarithmic order the Wilson coefficients and therefore 
BK(B — > X s u)l) are highly sensitive to the choice of the low-energy and matching scales 

3 More precisely, these are the expressions for the effective Wilson coefficients introduced in |lq]. 



A 



mainly due to dependence^ of (fi) on rj = a s (/j,w) / Oi s (/j,) . However, in the next-to-leading 
order the scale dependence of BR(£> — > Xgr^j) in general is significantly reduced |L8|, p(| [H], |22 



Cfl{^w) and Oj>%(p 



At the matching scale nw the coefficients Cj, i=l,2,3,5,6, are the same as in the Standard 

23|, [24fl. New Physics effects contribute only to C[ (/ijy), 



Model and can be found elsewhere EC 



0), 



w 1 



For Cfl{[iw) one has 



(0) v 



(o) H 
7,8 



|2^(0) ff 
I °7,8 



(2.2.4) 



(o) w (o) H ~ (o) H 

where Cf 8 ; , 8 y , Cf^ are given in ref . P4] . The superscripts W, H denote respectively the 
contribution from W-boson mediated and charged Higgs boson mediated loops. 
The direct CP-asymmetry in B — > X d ^ decay is defined to be 



a C p{B -> X d ~f) 



T{B^X dl )-T{B^X dl ) 
T(B^X dl ) + T(B^X dl y 



(2.2.5) 



It is known P5[-P5|, fl33| , fl35fl that the numerator of (2.2.5) starts at the order 0(a s ). Since 
the order 0{a 2 s ) corrections to (2.2.5) are known only partially, we expand (2.2.5) in a s keeping 
only the leading term. Consequently, we use one- loop RGE solution for a s {fi). 

The method of extracting acp(B — > X s ^ / y) in the SM and its extensions is described in 
ref's [I3|-[J5|, In the most general form the expression for B — > X d ^ decay asymmetry is 



(it is understood that the Wilson coefficients have to be taken at the low-energy scale /z) 



a C p{B 



-lm[Ci 0) *ci° h 



d \ d*l 

t 



20 Im[/ 27 ] 

— r H 

"1 2vr 




6tt 



27j C; (0) Im[C( 0) ]- 



+ 



Im[/ ; 

6tT 



27\ c '(0) 



Re(C. 



8 , 



(2.2.6) 



where r (the ratio of the absorptive parts of < d^i\0<i — \Oi\b > and < d r y\02 U — hOi u \b > 



matrix elements) is given in |35j and the bremsstrahlung coefficients fa and f'22 can be taken 



from [M\. The quantity can be easily derived from the expression for / 2 2 an d reads 



tic 
22 



(./: 



ctr 
22 



27 Jo 



r l/z 

/ dt(l - ztf 
Jo 



'G*(t) 




where z = m\ /ml, and G(t) can be found elsewhere [[3~5|, [17|, ^Ti 3. Note that our formula for 



the B — > Xd( s )7 rate asymmetry differs from that in ref. [45| by the last term in (2.2.6). This 
term arises from the interference of < d-yg\02 — \Oi\b > and < d-yg\02 U — \Oi u \b > matrix 
elements. 



4 c|°g(/^) depend (via C^°l(^w)) on the matching scale also due to dependence of the running top mass on 
the choice of fiw ■ However this dependence leads only to about 3% uncertainty in BR(£> — > X s / d \j) ]2(| . 



K 



The effective Hamiltonian for B — > X^( s )e + e decays can be written as 
H eff (b -> d( S )e + e~) = H eff (b -> d(a) 7 (+»)) 



H ' F Af s) (C 9 O 9 + C 10 O 10 ) 



where 



9 = d L (s L )7 M 6 L e7 M e, O 10 = d L {s L )^b L e^^e, 



(2.2.7) 



^o±cf (// 

,(0) 



CSV) 

(0) J 



ClO (/^w); 



(2.2.8) 



Cgm(pw) — C gw (flw) + |6l ^9,10 (A* 



and Pq and Cgi (a*w) are given in || ^9], |3(| |37fl . We consider the dilepton invariant mass 
in the range lGeV 2 < s < 6GeV 2 : then the long-distance effects arising from the contribution 
of cc and mi resonances can be neglected so that the decays B — ► X d ( s )e + e~ are well described 
by the transitions 6 — > d(s)e + e~ . The (partially integrated) direct CP-asymmetry in the decay 
B — > Xae + e~ is defined by 



a C p{B -> X d e + e" 



Afi( J B -> A d e+e~) - A£(£ -> X d e+e- 



(2.2.9) 



where A£(5 -> X d e+ e -) = / 7 2 " ^5(5 -> X d e + e~)/ds, s = s/m 2 b . The quantity dB[B -> 

A d e + e~)/<is can be expressed in terms of C7 (//) and Cg^(p) = Cg°^(/i)(l + w(s)) + 
where tu(s) and F 9 (s), q=s,d, are given in ref. 0. Denoting Cg '(/i) = 6 + [(A^A^*) / 1 2 ]^ 2 
one gets for i? — > X^e^e" decay rate asymmetry 



a C p 



ds(s 



-12 



M 2 
|A t | 2 



;i±2s)Im[£*£ 2 ]±6Re[Cf]Im[£ 



Im[C} 0) ]Im[6] - 12Im[C^ u; ]Im[e 



(0)1 



(2.2.10) 



J ds(s -!)*((! + 2s)[\C% ff \ 2 + \C w \ 2 ] ±12Re[Cf C e 9 ff *} + 4^|Cf | 2 



where the limits of the integrations are specified above. 



2.3 Constraints on parameters of theory. 



We use the following values of the parameters of the theory: M] 



m 



pole 



W 



80GeV, 



m 



pole 



(4.8 ± 0.15)GeV, mjm b = 0.29 ± 0.02 

0.0055 ±0.0013), 
The 



21, 24, B 



(174 ± 5)GeV, a s (M z ) = 0.119 ± 0.002 

(corresponding to r = 0.145i^f, Im[/ 27 ]/(27r) = 0.011 ±0.001, lm[f^f]/n = 
a em = 1/(130.3 ± 2.3) and BR sl = BR(B -> X c e + u) = (10.49 ± 0.46)% 
low-energy scale p and the matching scale pw are varied between rrib/2, 2m b and M w , m t 
respectively. 

In the Wolfenstein parameterization the necessary CKM-f actors A^,A? s) and V* are given 



by 



a; 



AA 3 (1 -p + ifj), 
-AA 2 (l-A 2 /2±A 2 (p -*??)), 



A\ 3 (p-ifi), 
AX\p- iV ), 



= Ay 



(2.3.1) 



where p = (1 - \)p, fj = (1 - ^77, A = sin0 c « 0.22 and A=0.819 ± 0.035. The restrictions 
on the parameters p and r\ can be obtained from the unitarity fits, which yield |48, |49| 



p 2 + v 2 = Q.423 ± 0.064, 38°<7<81° (2.3.2) 
16° < (3 < 34°, 75° < a < 121° (2.3.3) 

where a, j3 and 7 are the angles of the unitarity triangle. The conditions (2.3.2) lead to the 
following constraints on p and 77: 

0.05 < p < 0.38, 0.22 < 77 < 0.48 (2.3.4) 

During the numerical analysis we use (2.3.2) and (2.3.3) rather than (2.3.4) since it turns out 
that the correlations between the values of p and 77 reduce significantly the derived range for 
the CP-asymmetryQ. 

The restrictions on the parameters £ t and coming from the bounds on t- and b-quark 
Yukawa couplings have been discussed in subsection 2.1. It has been pointed out that while 
~ 1 or smaller, can be much larger than unity when tan/3 > 1 or tan/5 <C 1. This 
means that if no other constraints on the 2HDM parameters are imposed, the Wilson coefficients 
CfV) and Cf\p) can be significantly larger than in the Standard Model. The additional 
constraints on can be obtained from the experimental bound on B — > X s ^f 

branching. We require that 

BR 2HDM (< B -> X sl >) = BR exp (< B -> X sl >) = (3.15 ± 0.54) x 10" 4 (2.3.5) 

where BK 2HDM (< B—> X s 7 >) is calculated in the next-to-leading order according the formulae 
given in ref. pi |. Due to NLO effects the restrictions on \cf\p)\ coming from (2.3.5) can 



differ sizably from those given in [33]: for instance, in the case when sign (Re [C7 (//)]) = 
— sign(C , 7°' ) ' SA/ (yu)), \cf\p)\ can be up to 3 times larger than in the SM. Some deviation of 
\cf\p) \ from |C7°^ 5M (/i)| can also occur due to the fact that the SM prediction for the B —* X,7 
branching ratio only overlaps with the experimental data. 

As one can deduce from (2.2.4) and the relation |cf )H \ < |C} 0)H | |g§ |§, the condition 
(2.3.5) restricts the product £ t £b, whereas the parameter ^ and hence the Wilson coefficients 
Cg(ji) and C}q(ji) are still constrained only due to the bounds on t- and b- quark Yukawa 
couplings. Additional restrictions on (respectively on Cg°l (p)) can be obtained from the 
experimental measurements of the B — B mixing. For the central values of the parameters 
entering the expression for Am^ one obtains that < 0.3 when m#+ ~ lOOGeV |47| . However 
when taking into account experimental uncertainties of these parameters (especially important 
are those of Jb^Bb = (215 ± 40)MeV |[48|| ) one derives much weaker constraint. Thus, for 
m H + = lOOGeV, p = 0.19 and 77 = 0.37 we obtain that < 0.9. The restriction on \£ t \ again 
becomes stronger when p — > 0: |^ t | < 0.7 now. On the contrary, for p — > p max = 0.38 the bound 
on |£ t | coming from the experimental constraints on B — B mixing becomes inessential. 

The restrictions on the parameters £ t and ^ are also correlated with the charged Higgs mass. 
In our calculations we vary it as lOOGeV < < AOOGeV. For larger values of the 

5 Note that the 2HDM fits for the parameters p and r) can in general be different from those of the SM. This 
source of deviations from the SM predictions however is out of the scope of the present paper. 
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Table 3.1: a CP (B -> X dl ) in the Model I. 













a C p(B -> 


X d7 )/10" 2 










rriH+, GeV 






= M w 






= m t 








n = 


m b /2 


\l -- 


= m b 


= 2m b 


ji = m b /2 


V z 


= m b 


n = 


- 2m b 


100 


13 


+ 36 


11 


+ 33 


10^30 


13^32 


11 


+ 30 


10 


+ 27 


200 


13 


^35 


11 


^33 


10^30 


13^32 


11 


^30 


10 


+ 27 


400 


13 


-=-31 


11 


^28 


10^26 


13 -=-31 


11 


^28 


10 


+ 2Q 



theoretical predictions for BR(£> — ► X s ^) (and hence the restrictions on C7 (ji) and Cg (//)) 
become inaccurate, as the formulae for BR(I? — > Xg'j) contain terms proportional to log — — 



24 



H+ 



Some bounds on the 2HDM parameters could be obtained also from the experimental con- 
straints on the neutron electric dipole moment (NEDM). Note however that these constraints 
are known to suffer from the uncertainty of the hadronic matrix element < O g (/j, n ) > [ 501 , where 
Og = I f abc e Svalz 1 Ga/3^x8^Xu (here G denotes the gluon field strength tensor). When taking into 
account the uncertainty in the choice of < O g {ji n ) > (unlike (47| where a particular value of 
< Og(p, n ) > is used) the NEDM constraints restrict Im[£ t £&] weaker than the condition (2.3.5). 



3 Direct CP-asymmetry in the decay B —> X^. 

Though the condition (2.3.5) strongly restricts the 2HDM parameter space, the predictions 
of this model for the B — ► Xdj decay asymmetry can deviate significantly from the result 
derived in the Standard Model. This is mainly due to the fact that the condition (2.3.5) puts 
constraints on \Cj (//)| while its phase remains in general free. In particular, even being real 
cf\ix) and consequently acp{B —>■ Xd'j) can have the sign opposite to that in the SM. In other 
words the measurement of the direct CP-asymmetry in B — > Xd'j decay offers an alternative to 



the method proposed in [|5T], |31|] where the authors suggest to use B — > X s e + e branching and 
forward - backward asymmetry to determine the sign of Cf V)- The deviations of the 2HDM 
predictions from those of the SM can occur also due to the contribution of the Higgs doublets 
vacuum phase which gives rise the terms of (2.2.6) proportional to ImfCy ^/!)] and ImfCg '* (//)]. 
This contribution can lead both to the change of sign of a>cp and to the increasing or suppression 
of its absolute value as compared to that in the SM. So in the two-Higgs doublet extension of 
the Standard Model there are two sources of deviations of the predictions for acp(B — ► X^) 
from those of the SM. These sources are connected respectively with the change of sign of 
Re[C? (fi)], as compared to that in the SM, and with the presence of an additional (to the 
CKM one) weak phase in theory. In the Model III both of sources of deviations from the 
SM predictions are actual hence we investigate their interplay in details. On the contrary in 
the Models I and II both of these sources are irrelevant. Indeed, as it was discussed in the 
previous section, the parameters £ f , and hence C?(fi) and C$ \/j<) are real. Furthermore, 
for tan/3 > 1 the sign of Cf\n) is the same as in the SM (i.e. Cj°\jj) < 0). As a consequence 
in the Models I and II the predictions for a>cp(B — > X d -f) are close to those of the Standard 



Table 3.2: a CP (B -> X d 7) in the Model III. 
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Table 3.3: a C p{B -> X d7 ) in the Model III for m H + = lOOGeV, & = 0.5 and a) = 4 - the 
case when C7 (p) is real and positive; b) = 2±2i - the simplest case when arg(C? (fi)) — » 
±tt/2. 
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Model. 

Varying the parameters of the theory within the intervals specified in subsection 2.3 one 
finds that in the SM acp(B — > X^y) = (10 -4- 27)%. The uncertainty in the determination of 
acp is connected mainly with rather large dispersion of values of CKM parameters p and rj (for 
the detailed analysis of the behavior of acp(B — > Xd^y) with p and r] see ref. [031 ) ■ Indeed, if we 



fix these parameters e.g. at their "best fit" values, p = 0.19, rj = 0.37, then the uncertainty for 
the decay asymmetry is reduced significantly: one gets acp(B — > X d j) = (144 20)% now. The 
decay asymmetry shows moderate dependence (~ 15%) on the low-energy scale p. Surprisingly 
the dependence of a C p(B — > X^) on the matching scale pw is negligible. 

Note that obtained maximum (minimum) value of B — > X<fy decay asymmetry is about 1.3 
(1.5) times smaller (larger) than the one reported in ref. As it is mentioned above, this is 



mainly due to taking into account the correlations between the restrictions on p and rj. 

As it is noted above, in the Models I and II the contribution of the New Physics to the 
decay asymmetry is not sizable. In the Model I acp(B — > Xdj) = (10 4 36)% as one can see 
from Table 3.1. In the Model II one must take ra#+ > 400GeV when applying (2.3.5): then, 
as the numerical analysis shows, acp(B — > X^) remains within the SM interval. 

In the Model III the predictions for the B — > X^ decay asymmetry can differ from those 
of the SM significantly. As one can see from Table 3.2, acp(B — > X^) can be negative 
now. Furthermore it can be arbitrary small in absolute value and hence unmeasurable. Large 
deviations from the SM predictions concern principally the minimum value of B — > X^ decay 
asymmetry. As for the maximum value of CP-asymmetry, it can be larger than in the Standard 
Model only about 1.3 times. The deviation of a™p x (B — > A 7 ^) from the SM prediction is due 
to the fact that the experimentally allowed range for BK(B — > A s 7) only overlaps with the SM 
predictions. 



o 




Figure 1: Maximum and minimum values of acp{B — > X^y) as functions of 
MCf M)/Cf Ml for m H + = lOOGeV and fi = m b /2 (dashed line), \x = m b (solid line), 
\i = 2m b (dotted line). 



Thus we obtain that in the Model III a CP (B -> X d -f) = (-25 4- 37)%. Note that the 
maximum (absolute) value of CP-asymmetry is 1.5 - 2 times larger than in [[35|j . This is due 
to the use of new experimental results for the CKM parameters, taking into account the gluon 
bremsstrahlung effects and the scale uncertainty of the obtained results as well as due to the 
fact that our constraints on C\ (ji) differ from those of [35|| . 

We investigate in details the scale dependence of a>cp(B — > X d "f). As one can see from 
Table 3.2, minimum value of acp{B — > X^) is highly sensitive to the choice of low-energy 
and matching scales: it varies from —25% to —7% when varying /i and fiw between m b /2, 2m b 
and m t ,Mw respectively. In Table 3.3 we present the derived interval of B — > X,f) decay rate 
asymmetry for the fixed values of the parameters £ t and It is easy to see that large \i and 
dependence for a C p(B — > X d ^) is observed both when C?(fi) changes the sign, as compared 
to the SM, and when the effect of Higgs doublet vacuum phase becomes sizable. So, unlike the 
SM, in the 2HDM the direct CP-asymmetry in B — > decay can be highly sensitive to the 
choice of the renormalization scale when the New Physics contribution is sizable. 

Finally we examine the influence of the Higgs doublet vacuum phase on the obtained results. 
For this purpose we present in Fig. 1 acp(B — ► X d ^) as a function of lm.(C? (ji)/\C? 
Though the obtained results are renormalization scale dependent, they have the following com- 
mon feature. As one can see from Fig. 1, maximum (minimum) value of acp{B X^) 
derived for| lm(Cj°\fi)) = is close or coincides with the absolute maximum (minimum). On 

6 It is easy to see from (2.2.4) that in this case Im(Cg°' ) = as well. 

1 n 



the contrary, for lm.(C? ([i)/\C? ~ 1 the allowed interval of CP-asymmetry shrinks about 
1.5 times: a>cp{B — ► X^y) = (—17 -f- 20)% now. Thus one may conclude that in the 2HDM 
the new source of CP-violation either holds the allowed interval of B — > X^j decay asymmetry 
unaltered or leads to the suppression of \acp(B — > Xdj)\. 

It is important to stress that because of the Higgs doublet vacuum phase effects the B — > Xjy 
decay asymmetry can be arbitrary small in absolute value and hence unmeasurable. This is 
in contrast to the case with real Cf\p) and Cq°\p) when there is a lower bound on the 
B — > X rf 7 decay asymmetry: then acp(B — > X^) = (—24 -i — 3)% when being negative, 
and acp(B — > X^) = (10 4- 36)% when being positive. Notice also that the behavior of 
a>cp(B — > X^y) with Im^C? (/x)] is different from that of acp(B — > X s j). As it was shown in 



ref's. j38|, |35], p4[ , due to the presence of new source of CP- violation connected with spontaneous 
CP-breakdown the absolute value of acp(B — > X s ^) in the 2HDM can be a few times larger 
than in the SM. 

One can also expect that the behavior of acp{B — > X^n) with respect to any new source of 
CP-violation will be similar to that in the Model III for any extension of the Standard Model 
where the New Physics contribution to (p) is of the order or smaller than that to Cf^ (p) . 
Notice however that in the models with enhanced chromomagnetic moment operator the effect 
of new phases can be sizable ||45|| . Particular examples of such models are the supersymmetric 
theories. So far detailed analysis has been performed only for B —>■ X^ rate asymmetry 



56[ . The investigation of B — > X^ decay asymmetry in the MSSM has been carried out only 
for the case when the supersymmetry breaking parameters are real JIB] . It is interesting also 
to examine the case when these parameters are complex (and their phases are not suppressed). 
This work is in progress now. 



4 Direct CP- asymmetry in the decay B — > X^e + e . 

The deviations of acp(B — > X^e + e~) from the SM predictions can occur due to New Physics 
contribution to cf (/x), Cg° (/x), Cjo^/i). One must however take into account that the 2HDM 
effects modify these Wilson coefficients in different ways. One may use for C?V) the SM 
prediction since numerically Cg(fi) < 0. 025(7^ (fi). On the contrary the deviation of C^(p) 
from its SM value can be sizable. The Higgs-mediated loops contribution makes |Ci^(/i)| and 
hence the denominator of (2.2.10) larger than in the SM. Note however that the deviations 
of clo'(fi) from C$ SM {[i) depend on the restrictions on the parameter £ t and hence on the 
choice of the CKM parameters. These deviations are the largest for p — > p max and negligible 
for p — > Pmi n . As the numerical analysis shows, in the Standard Model the latter choice of 
p corresponds to the minimum value of B — > A^e + e~ rate asymmetry. As a result, the New 
Physics contribution to C{q (/V) holds the allowed interval of the acp(B — > X d e + e~) unaltered 
as compared to that in the SM (note that with the definition (2.2.10) acp{B Xde + e~) is 
positive in the SM). 

It is easy to see that if \Cr(n)\ ~ \Cf J (//)| < Cg (a*), the New Physics contribution to 
Cf\p) has negligible impact on B — ► X^e + e~ rate asymmetry. In the case when \cf\p)\ > 
1^7°^ (/-Ol (generally speaking in this case Re(Cf } (/i)) > 0) the denominator of (2.2.10) again 
increases. The numerator of (2.2.10) can be both larger and smaller than in the SM. Possible 
enhancement of the numerator cancels with the enhancement of the denominator - thus in the 
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Table 4.1: acp(B — > X d e + e ) in the SM for different choices of p and pw- 
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Table 4.2: a CP (B X d e + e~) in the Model III for m H + = (1004400)GeV and different choices 
of p and nw 
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2HDM the direct CP-asymmetry in B — > X^e + e~ decay in general does not exceed by more 
than (104 15)% its maximum value in the SM. However a>cp{B — > X^e + e~) can be significantly 
smaller than in the SM. Thus we disagree with the result of ref. [|37|, according to which the 
effects of New Physics in general make acp{B — > Xde + e~) larger than in the SM. 

Let us proceed to concrete results. The derived interval of acp(B — > X<ie + e~) in the 
Standard Model for different choices of the low-energy and matching scales is given in Table 4.1. 
As one can see, the direct CP-asymmetry in the decay B — > X^e + e~ is highly sensitive to the 
choice of the low energy scale: one gets ~ 65% uncertainty in the obtained results when 
varying /i from m b /2 to 2m b . This uncertainty can be much larger and exceed 100% when 
examining the renormalization scale dependence of acp{B — > X^e + e _ ) with the use of so- 
called Kagan-Neubert approach |26|| . It is interesting that varying the matching scale 
from My/ to m t , one obtains that the predictions for CP- asymmetry are shifted by the same 
magnitude independently of the choice of fi: by about (—0.7) x 10~ 2 for the minimum value 
and —(1.2 4 1.5) x 10~ 2 for the maximum value. 

For the fixed values of the low-energy and matching scales the uncertainty of acp{B — > 
Xde + e~) in the SM is connected mainly with large dispersion of the CKM parameters p and 
r], analogously to acp(B — > X d ^). For the "best fit" values p = 0.19, n = 0.37, a>cp(B — > 
X d e + e~) = (4.3 ±0.4)%. The behavior of B — > X d e + e~ decay asymmetry with p (equivalently, 
rj) for different choices of \J p 1 + n 2 and for p = m b , p w = M w is presented in Fig. 2. As one 
can see from this figure, a C p[B X d e + e~) increases as p or p' 2 + v 2 increases. 

The renormalization scale dependence makes the results for acp{B X d e + e~) in the SM 
ambiguous: when choosing p = 2m b and pw = M w , the direct CP-asymmetry in B — > X d e + e~ 
decay can reach 10%. On the other hand, for p = m b /2 and pw = m t , acp(B — > X d e + e~) is 
smaller than 1% thus being hardly measurable in future experiments. This ambiguity preserves 
also in the two-Higgs doublet extension of the SM. 

The predictions of the Models I and II coincide (within (10 4 15)% accuracy) with those of 
the Standard Model. For the Model III the derived interval of B — > X d e + e~ rate asymmetry 
for different choices of the low-energy and matching scales is presented in Table 4.2 . As one 
can see from this table, except of the case when p = m b /2 and pw = m t , the direct CP- 
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Figure 2: The dependence of a>cp{B — > X d e + e~) on CKM parameters for p 2 + i] 2 = 0.36 2 ; 0.42 2 ; 
0.48 2 (lines 1,2,3 respectively). Other relevant SM parameters are taken at their central values, 
the low-energy and matching scales are chosen as /i = m?,, nw = Mw- 



asymmetry in B — > X d e + e~ decay in the 2HDM does not exceed the SM maximum value more 
than (10 -r 15)%. On the contrary, the minimum value of acp{B — > X d e + e~) can deviate from 
the SM result significantly: it can be 1.5 -j- 2.5 times smaller than in the SM depending on the 
choice of the low-energy and matching scales. B — > X d e + e~ decay asymmetry can also change 
the sign, as compared to that in the SM. However in this case it remains small in absolute value 
and hence unmeasurable. 

Thus our main result in this section is that the 2HDM effects in general decrease the 
direct CP- asymmetry in B — » X d e + e~ decay as compared to the predictions of the SM. On 
the other hand, in the SM for the appropriate choice of the low-energy and matching scales, 
B — > X d e + e~ decay asymmetry can be driven to be arbitrary small (in absolute value). It 
remains small also in the 2HDM. As a result, the SM and 2HDM intervals of the CP-asymmetry 
do not differ significantly: one obtains acpiB — > X d e + e~) = (0.2 9.8)% in the SM and 
a CP (B -> X d e + e-) = (-0.5 -=- 10.7)% for the 2HDM. Thus when calculating B -> X d e + e~ 
decay asymmetry in the lowest non-vanishing order one is not able to conclude if this quantity 
is efficient for testing the physics beyond the Standard Model. To clarify this point one should 
include higher-order corrections to acp(B — > X d e + e~). 



5 Summary and Conclusions. 



We have considered the direct CP-asymmetry in the inclusive B — > Xjy and B — > X^e + e~ 
decays in the framework of general extension of the Standard Model with two Higgs doublets. 
The investigation has been performed in the lowest non-vanishing order in the renormalization 
group improved perturbation theory taking into account the restrictions on the parameters of 
the model coming from the assumption of the absence of New Physics between the electroweak 
breaking and unification scales as well as the experimental constraints from the measurements 
of the B — ► Xg'-f branching ratio and B — B mixing. 

It has been shown that in the most general version of the 2HDM (Model III) large deviations 
for acp(B — > X^) from its SM prediction are possible. In the case of the absence of a new 
CP- violating phase (additional to the CKM one) we have found that acp{B — > X^rf) can be 
both negative and positive varying in the ranges (—24^—3)% and (10-=- 36)% respectively. The 
measurement of the sign of acpiB — > X,fi) can determine the sign of the Wilson coefficient 
C?(fi). When an additional source of CP- violation is present our prediction for a C p(B — ► 
X d ^) is (—25 -T- 37)%. So the only new effect connected with the Higgs vacuum phase is that 
\acp{B —>■ Xdj)\ can be arbitrary small and hence unmeasurable. 

Our results suffer from the large dependence on the low-energy (// ~ m&) and matching 
(Hw ~ Mw, m t) scales when the contribution of the Higgs boson mediated loops is sizable. For 
instance, the minimum value of acp{B — > X^y) varies from —25% to —7% depending on the 
values of the scales. 

The problem of the renormalization scale dependence is actual also for B — > X^e + e~ decay 
asymmetry. Here the situation is more dramatic: though for some choices of the low-energy 
scale New Physics contribution decreases a>cp{B — > Xde + e~) significantly, its allowed interval 
does not differ essentially from that in the SM. In other words the calculations performed 
in the lowest non-vanishing order of the perturbation theory do not allow to conclude if the 
B — > Xde + e~ decay asymmetry is efficient for testing New Physics beyond the Standard Model. 
We expect that the problem of the renormalization scale dependence will be resolved after taking 
into account 0(a 2 s ) corrections to acp(B — > Xdj) and 0(a s ) corrections to acp(B — > Xde + e~). 
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